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Vn  important  optical  characteristic  of  seawater  is  the  ratio  of  the  light  scat-  /15' 
tering  index  a to  the  attenuation  index  e,  called  the  quantum  survival  parameter, 
which  in  many  cases  is  known  with  insufficient  confidence.  The  widely  usedmgthod 
of  determination  of  this  quantity  involving  separate  measurement^  of (Qj 
frequently  associated  with  large  instrumental  errors^  Of  interest  in  this  connec- 
tion is  the  possibility  of  calculating  the  parameter  nvfrom  light  field  character- 
istics measured  directly  in  a scattering  medium.  ' 


Described  below  is  a method  of  determining  th^  parameter  A (a)  in  a plane-stra- 
tified inhomogeneous  medium  from  the  parameters  of  the  light  field  inside  the  medium 
when  its  surface  is  uniformly  illuminated  by  a flux  of  solar  radiation,  and  (b)  in  a 
homogeneous  medium  from  the  light  field  parameters  of  an  isotropic  source  placed 
inside  the  medium  at  a sufficient  distance  from  its  boundaries. 


In  describing  the  process  of  propagation  of  solar  radiation  in  the  sea,  it  is 
convenient  to  use  the  radiation  transfer  equation  in  the  following  form:^ 


cosO  ?)-!-»  Jz(«.  o'.  ?-?’)"  (O'. 


(i) 


where  B(0,<J>)  is  the  luminance  of  the  radiation  in  the  direction  characterized  by  /158 

angle  0 with  the  vertical  and  azimuth  <{>;  ~ 4)*)  is  the  indicatrix  of  single 

scattering  of  a ray  following  the  direction  (0,,<J>')  and  falling  on  the  scattering  volume 
in  the  direction  (0,$);  dui'  = sin  0 * d0 * d<|) * ; the  z axis  points  vertically  downward. 


In  a plane-stratified  medium,  the  luminance  remains  unchanged  in  the  direction 
0 = ir/2.  Integrating  Eq.  (2)  with  respect  to  <J>  and  choosing  the  corresponding  order 
of  integration  in  the  triple  integral,  we  obtain 


\4r]x 
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notation 
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It  follows  from  F,q.  (3)  that  15^(0)  is  the  azimuth-averaged  luminance  value  pro- 
pagated at  a given  depth  at  angle  0 with  the  vertical.  The.  quantity  x<j,(0)  is  the 
value  of  the  scattering  indicatrix  of  a ray  falling  on  the  scattering  volume  at 


Numbers  in  the  right  margin  indicate  pagination  in  the  original  text. 


angle  0 with  the  vertical,  averaged  over  all  the  directions  lying  in  the  horizontal 
plane.  X<J>(9)  is  independent  of  <J>'»  since  a periodic  function  of  (<p  - <p’)  is  integrated 
in  Eq.  with  respect  to  <f>  over  the  segment  [0,2tt].  The  initial  function  for  cal- 
culating X(j>(9)  is  the  experimentally  determinable  indicatrix  of  single  scattering  x(Y)- 

Taking  Eqs.  (3)  and  (A)  into  account,  we  represent  (2)  in  the  form 


where 
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Formula  (5)  makes  it  possible  to  calculate  A in  a plane-stratified  medium  at 
any  depth,  since  there  are  no  restrictions  on  the  shape  of  the  luminance  solid. 

In  the  case  of  symmetry  of  the  light  field  with  respect  to  the  vertical  axis 
(the  luminance  is  independent  of  the  azimuth),  expression  (6)  simplifies  to 


a 
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This  condition  is  fulfilled  in  the  sea  in  cloudy  weather  at  any  depth  or  in 
sunny  weather  when  the  sun  is  at  the  zenith,  and  also  in  the  state  of  optical  equi- 
librium under  any  conditions  of  surface  illumination.  The  condition  of  light  field 
symmetry  does  not  necessarily  imply  the  advent  of  the  state  of  optical  equilibrium; 
the  latter  may  fail  to  take  place  in  a vertically  inhomogeneous  medium. 

In  deriving  the  formula  for  calculating  A from  the  parameters  of  the  light  field 
of  an  isotropic  source,  we  will  use  the  transfer  equation  in  spherical  coordinates 

r)+.J  /.(#,  o’.  ?-?’)« (O’.  »■)-£-.  (7) 


where  r is  the  distance  from  the  source;  0 is  the  angle  between  the  specified  direc- 
tion and  the  direction  of  the  source;  <f>  is  the  azimuth.  The  luminance  B in  the  field 
of  the  isotropic  source  is  independent  of  the  azimuth  (the  luminance  solids  possess 
axial  symmetry  with  respect  to  the  direction  of  the  source).  We  will  consider  the 
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change  in  luminance  in  the  plane  perpendicular  to  the  direction  of  the  source:  0 = -j. 

The  left-hand  side  of  Eq.  (7)  may  be  represented  as 

'IB  (».  r > _ Ja  . 'I'  i **  . A.  /*) 
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It  follows  from  geometric  relationships  that  d&  = -rdO;  dr  ~ 
the  second-order  term,  we  obtain  from  Eq.  (8) 


-d?.d0.  Omitting 
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and  Eq.  (7)  becomes 
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Dividing  both  sides  of  the  equation  by  B(y)  and  using  the  notation  of  (4)  and  (6a),  /16 

and  considering  that  e = K + a,  where  < is  the  absorption  index,  we  arrive  at  the 
following  expression  for  A: 


or,  in  view  of  the  notation 
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When  r -*■  00 , Eq.  (13)  becomes  Eq.  (5). 

Finally,  when  0=0,  it  follows  from  Eq.  (7)  that 
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Also  using  the  expression-1  for  K 
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is  the  mean  cosine  of  the  radiation; 
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we  obtain  from  (14)  the  following  formula  for  calculating  A: 
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In  the  state  of  optical  equilibrium  «B(0)  3 ®F  = a»  and  E<I-  (21)  becomes 
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Thus,  four  formulas  have  been  proposed  for  calculating  A:  in  a plane-stratified 
medium  Illuminated  by  a wide  light  beam,  at  any  depth,  (5);  in  the  state  of  optical 
eqiulibrium,  (22);  from  the  light  field  of  an  isotropic  source,  (13)  and  (21). 

2.  Measurements  of  the  luminance  B(6,r)  of  the  light  field  of  an  isotropic  source 
and  luminance  of  the  light  field  of  solar  radiation  in  the  state  of  optical  equilibrium 
were  performed  in  the  Black  Sea  (1970),  Mediterranean  Sea  (1966),  and  certain  regions 
of  the  Atlantic  and  Pacific  Oceans  during  the  fifth  cruise  of  the  research  ship 
DMITRIY  MENDELEYEV  (1971).  The  instruments  and  methods  used  for  measuring  the  light 
field  characteristics  of  an  isotropic  radiator  are  described  in  Ref.  3.  "Deep" 
luminance  solids  were  measured  in  the  Black  Sea  by  using  the  same  instruments,  and 
in  the  other  waters,  by  means  of  an  underwater  photoelectric  luminance  meter  designed 
by  A.  S.  Suslyayev,  mounted  on  an  underwater  photometric  bench  with  a scanner. 

The  spectral  range  of  the  luminance  measurements  was  495±35  nm.  The  volume 
scattering  function  B(y)  was  measured  in  relative  units  by  a standard  method.  Data 
on  the  measurements  of  B(0,r)  and  B(Y)  in  the  Black  Sea  are  given  in  Tables  1 and  2. 
Table  1 gives  the  angular  luminance  distributions  in  the  light  field  of  an  isotropic 
source  at  different  distances  from  the  latter,  and  also  in  the  state  of  optical  equi- 
librium. The  values  of  log  — are  given,  where  is  the  light  flux  of  the  source 

Table  2 gives  values  of  the  volume  scattering  function  B(y)  ~ not  normalized,  in 
relative  units  - from  measurements  at  st.  100  (Black  Sea).  Calculations  of  the 
quantum  survival  parameter  A from  formulas  (5),  (13),  (21),  and  (22)  were  made  on  /16 
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a Minsk-22  computer.  The  results  of  the  calculations  of  A for  different  areas  of 
the  World  Ocean  are  summarized  in  Table  3,  which  shows  that  the  four  formulas  for 
calculating  A at  st.  100  (Black  Sea,  1970),  where  the  most  complete  measurements  of 
the  light  fields  were  performed,  gave  very  consistent  results.  The  discrepancies  in 
the  value  of  A do  not  exceed  the  limits  of  error  estimated  below.  The  much  greater  /16 
discrepancies  observed  at  st.  367  and  397  may  be  explained  by  the  following  factors: 

(a)  measurements  of  "deep"  luminance  solids  were  made  during  the  day,  whereas  the 
distribution  of  luminance  in  the  light  field  of  an  isotropic  source  was  measured  at 
night  (the  ship  could  have  been  drifting  into  waters  with  somewhat  different  optical 
characteristics);  (b)  the  values  of  the  parameter  A,  calculated  from  formula  (13), 
were  averaged  over  a certain  depth  range,  whereas  formulas  (5)  and  (22)  determine  A 
a£  a depth  of  100  m.  The  highest  value,  A = 0.82,  corresponds  to  the  Black  Sea,  and 
the  lowest,  to  clear  waters  of  the  Sargasso  Sea  in  its  southern  fringe. 


It  is  interesting  to  compare  the  value  A = 0.32  obtained  for  the  Sargasso  Sea 
(at  a depth  of  300  m)  with  the  value  of  this  parameter  for  "ocean  waters  of  maximum 
purity,"  proposed  in  Ref.  4.  According  to  the  data  of  this  study,  in  ocean  water 
of  maximum  purity  A = 0.32  at  X = 550  nm. 


The  last  two  columns  of  Table  3 give  values  of  the  attenuation  index  e and 
mean  scattering  angle 

f iMU** 
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characterizing  the  elongation  of  the  single  scattering  indicatrix. 

The  three  quantities  - quantum  survival  parameter  A,  attenuation  index  e and 
scattering  indicatrix  - which  completely  determine  the  intrinsic  optical  charac- 
teristics of  a medium, are  mutually  independent.  However,  for  natural  waters,  not 
all  combinations  of  (A,  e,  y)  are  equally  probable,  and  a relationship  is  observed  /165 
between  them,  as  illustrated  in  the  figure. 

The  figure  gives  values  of  A,  £ and  y from  the  data  of  our  measurements.  A 
smooth  decrease  in  A and  increase  in  y with  decreasing  index  e is  observed.  The 
following  explanation  of  this  dependence  is  possible.  On  passing  from  polluted 
waters  to  purer  ones,  primarily  the  concentration  of  the  coarsely  disperse^  suspen- 
sion decreases,  so  that  the  scattering  indicatrix  becomes  less  elongated  (y  increases), 





Table  3 


Region  of  measure- 
ment, No.  of  sta- 
tion, range  of 
J — ths 


Value  of  A averaged  over  the  indi- 
cated depth  range 


Calculation  from 
light  field  of 
isotropic  source 


Calculation  from 
light  field  of 
solar  radiator 


Black  Sea  st.  100, 
30-150  m 

Black  Sea,  st.  33, 
100-150  m 

Pacific  Ocean, 
Equator,  Cromwell 
Current,  st.  367, 
30-60  m 

Mediterranean  Sea, 
st.  34,  100-150  m 

Pacific  Ocean,  re- 
gion southeast  of 
Japan,  st.  410, 
30-150  m 

Pacific  Ocean,  re- 
gion of  high  trans- 
parency near  Kuka 
Islands,  st.  397, 
30-150  m 

Pacific  Ocean, 
Equator,  Cromwell 
Current,  st.  367, 
120-200  m 

Atlantic  Ocean, 
southern  Sargasso 
Sea,  st.  347, 

30-100  m 

Atlantic  Ocean, 
southern  Sargasso 
Sea,  st.  347, 
200-300  m 


Formula  Formula  Formula  Formula 

(13)  (21)  (5)  (22) 
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Relationship  between  A,  e and  Y based  on  results  of  measurements  in  different 
water  areas  of  the  World  Ocean  (see  r.lso  Table  3). 

Arabic  numerals  next  to  points  ^Y,  deg;  1 - optical  properties  of  pure  water 
(Rayleigh  scattering,  absorption  only  by  water  molecules). 


In  the  three-dimensional  space  formed  by  the  quantities  A,  e and  y,  the  optical 
properties  of  the  investigated  waters  fall  near  some  curve  (see  figure)  apparently 
corresponding  to  the  zone  of  combinations  of  these  quantities  most  frequently  en- 
countered in  natural  waters.  Hence,  from  one  of  the  three  indicated  characteristics 
one  can  estimate  the  other  two.  This  hypothesis  opens  up  the  possibility  of  classi- 
fying the  waters  of_the  World  Ocean  according  to  the  location  of  the  point  with 
coordinates  (A,  e,  y)  on  the  indicated  curve.  The  transition  from  the  classification 
of  waters  into  groups  differing  sharply  in  optical  properties,  as  was  done  in  Ref.  5, 
to  a continuous  distribution  of  waters  according  to  (A,  e,  y)  along  the  indicated 
curve  is  more  convenient  in  many  respects.  This  hypothesis  should  be  checked  by 
using  more  extensive  experimental  material  than  what  is  available  to  us  at  the  present 
time,  and  the  scatter  of  experimental  values  about  the  specified  regression  line  should 
also  be  quantitatively  estimated. 

3.  The  error  in  the  determination  of  A from  formulas  (5),  (13),  (21)  and  (22)  / 1 66 

is  due  to  experimental  errors  in  the  measurements  of  luminance  and  volume  scattering 
function,  which  were  used  as  the  basis  of  the  calculation,  and  also  to  errors  of  the 
numerical  integration.  In  the  calculation  with  the  Minsk-22  computer,  the  integration 
step  chosen  was  small  enough  so  that  the  computational  error  did  not  exceed  a fraction 
of  a percent;  it  may  therefore  be  neglected  in  comparison  with  the  measuring  error. 

The  error  in  the  luminance  measurement  was  in  the  range  of  ± 18%  for  0°  ^ 0 < 10°;  12% 
for  10°  < 0 $ 20°;  9%  for  20°  < 0 $ 30°;  7%  for  30°  < 0 < 60°;  6%  for  60°  < 9 « 180°. 

In  the  measurement  of  3(y).  the  error  did  not  exceed  ±1.0%  in  the  angular  range  0°  < y 
10°  and  ±7%  at  y > 10°. \ 

The  limiting  error  ln\h  was  estimated  from  a formula  derived  from  F!q.  (13): 
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with  ~ = 0.1,  as  was  shown  in  Ref.  3;  Ar  = 0.5  m was  the  scale  division  of  Lhe  winch 
K Aon, 

Indicator;  = 0.1  is  the  maximum  deviation  from  the  mean,  obtained  by  processing 

“b 

a large  number  of  measured  functions  B(0);  -r  = 0.06  is  the  maximum  deviation  in  I, 


calculated  as  follows.  In  the  calculation  of  I from  Eq.  (12),  the  absolute  values 
of  the  functions  B(0)  and  8(Y)  do  not  affect  the  result  - only  their  relative  variation, 
i.e.,  the  "degree  of  elongation,"  is  important.  Having  specified  the  above  limiting 
errors  in  B(0)  and  8(Y)  * we  drew  the  curves  B'(0),  B"(0)  and  8' (y) » 8"(y)  within  the 
confines  of  the  error  zone  in  such  a way  that  B'(0)  and  8'(y)  would  have  the  maximum 
elongation,  and  B"(0)  and  8"(y)»  the  minimum  elongation  along  the  axis  of  ordinates. 
Calculation  of  I values  for  any  two  pairs  of  the  indicated  functions  gives  four  values 
of  the  integral  I.  This  procedure  was  carried  out  for  three  distances  r from  the 
source  by  using  the  data  of  measurements  at  st.  100  (in  the  Black  Sea).  The  maximum 
deviation  of  the  integral  calculated  in  this  manner  for  the  entire  set  of  calculated 

cases  was  6Z  of  the  value  of  the  integral  itself;  therefore,  ~ = 0.06  was  substituted 

into  Eq.  (23)  as  the  limiting  value. 

•i 

Also  substituted  into  formula  (23)  were  values  of  the  characteristics  corres- 
ponding to  measurements  in  the  Black  Sea:  a = 0.9,  < = 0.03  m~l,  r = 100  m,  I = 1.28;  /If 

as  a result,  AA  = 0.065  was  obtained.  According  to  measurements  in  the  Sargasso  Sea, 

a = 0.77,  ic  =*  0.014  r = 100  m,  I = 4.47;  correspondingly,  = 0.062.  Considering 

that  these  cases  are  extreme,  the  errors  in  A values  calculated  from  the  data  of  the 
remaining  measurements  may  be  estimated  as  being  approximately  the  same.  It  is  inter- 
esting to  note  that  the  errors  in  <,  r and  a make  a very  small  contribution  to  the 
error  of  A,  owing  to  the  fact  that  these  quantities  are  present  in  both  the  numerator 
and  the  denominator  of  formula  (13).  The  determining  error  is  the  one  in  I,  which 
we  estimated  in  the  manner  described  above.  On  the  whole,  the  limiting  error  of  the 
above  A values  in  the  calculation  using  Eq.  (13)  may  be  considered  equal  to  ±7%  of  A. 

The  corresponding  estimates  were  made  for  the  calculation  using  formulas  (5),  (21) 
and  (22);  the  limiting  errors  given  in  Table  3 incorporate  these  estimates. 

It  is  our  view  that  the  proposed  method  gives  a more  accurate  value  of  the  para- 
meter A than  the  other  methods  in  use  at  the  present  time,  since  the  calculation  in- 
volves the  use  of  data  of  photometric  measurements,  which  do  not  require  calibration 
of  the  instruments  (brightness  and  indicatrix  meters)  in  absolute  units.  The  dimen- 
sionless coefficient  A is  determined  on  the  basis  of  measurements  of  the  two  dimen- 
sionless functions  b(0)  and  8(Y) * this  is  preferable  to  determining  A by  dividing 
by  each  other  the  two  dimensional  quantities  o and  e,  measured  independently  in 
absolute  units  with  various  instruments  and  by  different  measuring  techniques. 

Concl us  Ions 

1 . A method  is  proposed  for  determing  the  light  quantum  survival  parameter  A in 
seawater  from  the  characteristics  of  the  light  field  of  an  isotropic  source  and  light 
f icld  of  radiation  in  the  sea.  This  parameter  was  calculated  for  several  .‘trcas  of  the 
World  Ocean.  The  highest  value  A = 0.8210.06  for  the  waters  studied  was  observed 
la  the  Black  Sea,  and  the  lowest,  0.3210.03,  in  the  Sargasso  Sea. 

A relationship  was  observed  between  the  three  intrinsic  characteristics  of  sea- 
waters - the  parameter  A,  attenuation  index  £,  and  shape  of  the  scattering  indicatrix, 
characterized  by  the  scattering  angle  Y*  The  remaining  characteristics  of  the  waters 
studied  fell  on  a single  curve  in  the  space  of  (A,  e,  >),  which  opens  up  new  possi- 
bilities for  classifying  seawaters. 

The  authors  are  grateful  to  0.  V.  Kopelevich  and  V.  M.  Pavlov  for  kindly  supplying 
the  data  of  measurements  of  the  volume  scattering  function  3(y)  tinting  the  fifth 
cruise  of  the  research  ship  DMITRIY  MENDELEYEV. 
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